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I. INTRODUCTION 

This work completes the formalism that was partially developed in Ref. [ 1 ] . In that work we determined the finite- 
volume spectrum of a relativistic scalar quantum field theory by calculating an appropriate finite-volume correlator 
to all orders in perturbation theory. All contributing Feynman diagrams were decomposed into products of finite- and 
infinite-volume quantities. Summing these factored expressions gave a result for the correlator in terms of finite-volume 
kinematic functions and infinite-volume scattering quantities. The poles in the resulting expression determine the 
energies of the finite-volume states, and identifying their locations gave a relation between the finite-volume spectrum 
and infinite-volume scattering quantities. 

The central drawback of the result of Ref. [T] is that it depends on a nonstandard infinite-volume three-particle 
quantity, a modified three-particle K-matrix denoted /Cdf,3- For our formalism to be useful, it is necessary to relate 
^df,3 to the physical three-particle scattering amplitude AI3. This is what we achieve in the present work. This 
shows that, in a relativistic context, the finite-volume spectrum of three particles is determined by infinite-volume 
observables (up to corrections falling faster than any power of 1 /T, where L is the linear extent of the finite volume). 
This result was previously established using nonrelativistic effective field theory [ 2 ]. 

To give a sense of what is involved in relating /Cdf,3 to AI3, we recall from Ref. [T] the major differences between 
the two quantities. Both are obtained by summing contributions from diagrams with six external lines. However, 
divergences that are present in Mz are removed from JCdf,3 (leading to the subscript “df” for “divergence free”). 
These divergences are due to diagrams containing pairwise scatterings, such as that shown in Fig. [l] They have 
nothing to do with bound states but are instead a result of on-shell three-particle intermediate states I3HS]. It is not 
surprising that the finite-volume spectrum is related to a modihed quantity in which these divergences are removed. 
We stress that the terms removed depend only on the on-shell two-to-two scattering amplitude, A 42 - These can thus 
be separately computed and added to the divergence-free quantity to recover AI3. 



FIG. 1. Example of a singular contribution to the on-shell three-to-three scattering amplitude. Dashed lines are on-shell, 
amputated, external propagators, while the solid line is a fully dressed propagator. Filled circles represent two-to-two scattering 
amplitudes. The momentum flowing along the internal (solid) line goes on shell (leading to a divergent scattering amplitude) 
for particular choices of external momenta corresponding to two isolated two-to-two scattering events. 


A second and more important difference between AI3 and /Cdf,3 concerns the pole prescription used to define 
momentum integrals over products of propagators. In AI3 the standard ie prescription is used, while Aildf,3 uses a 
modified principal value prescription (denoted PV). The use of a nonstandard pole prescription is required in order 
to remove the unitary threshold cusp in the two-particle scattering amplitude, as was pointed out in Ref. [ 5 ]. As 
described in detail in Ref. [ 1 ] and sketched in Fig. such cusps, if not removed, generate important finite-volume 
effects when two-to-two scattering is considered as a subprocess of three-to-three scattering. Our prescription avoids 
such cusps, but at the cost of introducing a nonstandard infinite-volume quantity. 

A third difference betw e en Ad 3 and /Cdf,3 is that the latter is defined using an ultraviolet cut-off (the function 
H{k) defined in Eqs. (A 8 |-(AI 0 ) of the appendix) while AI3, being a physical quantity, is independent of the cut-off 
function. 

The final difference is that /Cdf,3 involves contributions in addition to the standard Feynman diagrams, which we 
denote “decorations”. These make the definition of /Cdf,3 very complicated, as explained in Ref. [T]. 

In this work, we relate these two quantities using an indirect approach. By extending the methods of Ref. [ 1 ], 
we obtain a relation between which is a finite-volume version of AI3, and /Cdf,3j Eq. (68) below. The desired 

relation between AI3 and /Cdf,3 is then reached by taking the infinite-volume limit. This approach avoids the need 
to directly reference the very complicated definition of /Cdf,3- As a side benefit, we obtain some intuition about the 
properties of /Cdf,3 itself. 

Combining the results of Ref. [T] with the present article gives a complete relation between the finite-volume energy 
spectrum and two-to-two and three-to-three scattering amplitudes. The result relies on no assumptions concerning 
the details of the underlying theory, no specific effective field theory and no power-counting scheme. Further, the 
relation is derived in fully relativistic held theory, and thus does not rely on a small momentum expansion. However, 
as it stands, the formalism is limited in two important ways. First, it can only describe systems of identical scalar 
particles with a Z2 symmetry that prevents the coupling of states containing even and odd numbers of particles. In 
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FIG. 2 . (a) Schematic plot of the the s-wave two-to-two scattering amplitude, Al2,s, as a function of two-particle CM energy, E2, 
showing the well known cusps at threshold {E2 ~ 2 m) in both real and imaginary parts, (b) Example of two-to-two scattering 
as a subprocess of three-to-three scattering. Because total energy and momentum are fixed, the energy and momentum of the 
top two propagators varies with the momentum of the bottom propagator. Thus the sum over k has important finite-volume 
effects arising from the two-particle unitary cusp, (c) In contrast to M2, IC2 is a smooth function at threshold. Working with 
K.2 in place of M2 resolves the issue of finite-volume effects from threshold cusps. 


particular, no two-to-three scattering is included. Second, there can be no two-particle resonances in the energy range 
considered. If the center of mass (CM) energy of the three-particle system is E*, then resonances must not appear 
below two-particle CM energies of E* — mQ We make this restriction more precise in the paragraph following Eq. Q 
below. There is an additional, less restrictive condition, namely that E* must lie above the single-particle pole and 
below the hve-particle production threshold (m < E* < 5 m). This is analogous to the restriction of the two-particle 
quantization condition to energies below the inelastic threshold. 

Despite these restrictions, there are several physical systems to which our formalism can be directly applied. In the 
context of QCD, one must consider three-meson systems with a quantum number that enforces the Z2 symmetry. For 
example, in the limit that isospin is an exact symmetry, states with even and odd numbers of pions decouple due to 
G parity. In this case one can project onto states of definite / and I3 and effectively obtain a system of three identical 
scalars. A simple example is the system with I = I3 = 3 . For all choices of isospin, our formalism is directly 

applicable up to the five pion threshold. The resonances that are present in some channels, e.g. the p if there is an 
7 = 1 two-particle subchannel, occur at energies above the five pion threshold]^ 

Another application of our present formalism is to study three kaon and three D meson systems, for which the 
conservation of strangeness and charm in QCD prevents coupling to other states. 

Finally, we note that our formalism also applies for constituent particles such as ultra-cold atoms and molecules. 
Here a non-relativistic description is usually expected to work very well so that the relativistic machinery we employ 
would appear to be overkill. However, in certain cases the relativistic result can give interesting information even 
very close to threshold. For example in a recent work. Ref. [ 7 ], we have studied the threshold shift of two- and 
three-particle states in \(j>^ theory. We expand the shift from threshold in powers of inverse box length, i.e. powers 
of 1 /L, and find that at 1 /L® the result differs from previous non-relativistic analysis [ 5 ]. This is also the order at 
which three-body couplings first arise in the expansion. Thus, in this context, three-particle and relativistic effects 
have the same volume scaling and it is natural to include them together. 

Going beyond perturbative expansions, we also note that our formalism allows for arbitrarily large two-particle 
scattering lengths. In particular, our results can be used to describe systems at unitarity, where the two-particle 
scattering length diverges]^ In this vein an interesting study has recently been published in Ref. [ 3 ]. Here the authors 
consider a non-relativistic three-particle system with two-particle scattering at unitarity. They then analyze how a 
shallow three-particle bound state would be shifted due to finite-volume effects. Their result is a clear quantitative 
prediction for the leading finite-volume shift, and reproducing this with the present general formalism would be an 
interesting and important check. 

The remainder of this article is organized as follows. In the following section, we recall the results of Ref. [ 1 ] that are 
needed here. In Secs. HI and IV we derive the relationship between and /Cdf,3. The former section includes only 

two-particle Bethe-Salpeter kernels, while the latter adds in three-particle kernels. In Sec. we take an appropriate 


^ Since E* is the CM energy for three particles, E* — m is the maximum possible energy for a two-particle subsystem. 

^ We caution that for sufficiently heavy pions the p resonance is pushed below the four particle threshold of the two-to-two system, and 
thus below the five-particle threshold in the three-to-three system. In this case the rho mass, rather than the 5m threshold, becomes 
the upper kinematic bound. 

® The scattering length a and effective range r are defined via 

p cot 5 (p) = -^— rp^ + O (p^), 

a 2 

where S is the s-wave scattering phase shift and p is the magnitude of three momentum for one particle in the two-particle CM frame. 
Our formalism is only valid if p cot 5{p) does not pass through zero for a particular p-value in the energy range considered. In the context 
of the unitary limit, this requires setting r and all higher order coefficients to zero before sending a —> 00 . 
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infinite-volume limit of M.z^l and obtain the desired relation between and /Cdf.a- We then show, in Sec. VI 
how to invert this relationship so that /Cdf,3 can be obtained from Section VII shows how the general relation 
between these quantities simplifies in two approximations studied in Ref. [T]. We conclude in Sec. VIII We include 
an appendix which collects the definitions of key quantities from Ref. [T] . 


II. FINITE-VOLUME CORRELATOR FROM PREVIOUS WORK 

In this section we summarize the relevant results from Ref. [T] , and rewrite them in a form useful for our subsequent 
analysis. We also introduce the finite-volume three-to-three scattering amplitude the quantity that plays a 

central role in our analysis. 

We begin by reviewing the set-up for the calculation presented in Ref. [T]. We work with a quantum field theory 
containing a real scalar, (j), having physical mass m. The theory is assumed to have a Z2 symmetry so that only even 
powers of (j) appear in the Lagrangian and thus even- and odd-particle states have zero overlap. The derivation of 
Ref. [T] also requires that the total center-of-mass (CM) energy is such that only three-particle states can go on-shell, 
and also such that there are no resonances in the two particle subchannels of the three-particle states. The theory is 
otherwise assumed to be completely general. In particular, it need not be renormalizable and can contain all possible 
even-legged vertices with no assumptions about relative coupling strengths. 

This theory is then considered in a finite, cubic spatial volume of extent L with periodic boundary conditions. The 
time extent is taken infinite. In this geometry, only spatial momenta whose components are integer multiples of 27 r/L 
can propagate. The volume extent L is also assumed to be large enough that exponentially suppressed corrections 
of the form can be neglected. Such terms are dropped throughout the derivation. Power-law corrections 

proportional to powers of lj{mL) are kept to all orders. 

The analysis of Ref. [T] considers a finite-volume correlator defined by 

Cl{E,P) = J d^xP^^^°-^-^\0\Taix)a\0)\0) . (1) 

Here T indicates time-ordering while cr{x) and cr^(a;) are interpolating fields coupling to states with an odd number of 
particles. The Fourier transform restricts the intermediate states contributing to the correlator to have total energy 
E and momentum P £ { 2 ttIL)Z^. Thus the energy in the center of mass (CM) frame is E* = Ve'^ — P'^. Our 
analysis concerns the region m < E* < 5 m, where only three-particle states contribute. Note that, although we 
work in Minkowski space, the finite-volume spectrum is not affected by this choice. Our result is applicable even 
if the energies are determined in some other way, for example by using a Euclidean correlator as in a lattice QCD 
calculation. 

The finite-volume correlator, Cl, has the skeleton expansion shown in Fig. This is explained in Ref. [T] (where 
precise definitions of the kernels B2 and B3 can also be found). Here we only note the essential features. Note 
first that Cl differs from its infinite-volume counterpart. Coo, only in that spatial loop momenta are summed rather 
than integrated. The key observation is then that replacing sums with integrals leads to power-law finite-volume 
corrections only when all particles in a cut can go on shell. Since we have chosen m < E* < 5 m, this can occur only 
for three-particle cuts. For loops with cuts involving five or more particles one makes only exponentially suppressed 
errors when replacing sums with integrals. As we neglect such corrections throughout, the Bethe-Salpeter kernels (for 
which cuts contain either one or else five or more particles) can be replaced by their infinite-volume counterparts. 
The same holds true for the dressing functions of the propagators. Thus the only impact of working in finite volume 
is that the propagators contained within the dashed rectangles in the figure have spatial momenta that are summed. 

The task accomplished in Ref. [T] is to factorize the power-law finite-volume dependence contained in Cl from 
quantities defined in infinite volume. The end result is 

CL{E,P) = CooiE,P) + iA' J: F,A. ( 2 ) 

I -I- Fa/Vdf.s 

Here A', A and /Cdf,3 (as well as Coo) are infinite-volume quantities, while F3 depends on L [as well as on E, P and 
the two-particle scattering amplitude—see below]. The result is a matrix equation, in which A' is a row vector, 
Adf,3 and F3 are matrices and A is a column vector. The index space is a direct product of two-particle angular 
momentum, parametrized by spherical harmonic indices £ and m, together with a discrete finite-volume momentum, 
k £ { 2 t: jL)'IA. So, for example, /Cdf,3 = ICd{,3-,k'e'm'\kim where k' and k are shorthand for discretized three-vectors. 

As we will see below, Coo, A and A! do not appear in the main results of this paper and so we do not recapitulate 
their definitions from Ref. [T]. As for /Cdf,3, the main task of this paper, as noted in the introduction, is to come 
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FIG. 3 . Skeleton expansion for the finite-volume correlator. Outermost blobs in all diagrams represent functions of momentum 
that are determined by the interpolating operators a and cr^. Insertions between these functions having four legs represent 
two-to-two Bethe-Salpeter kernels, while insertions with six legs represent the analogous three-to-three kernels, iB^. Lines 
connecting kernels and cr-functions represent fully-dressed propagators. As explained in the text, the kernels and dressed 
propagators can be replaced by their infinite-volume counterparts (in which internal loop momenta are integrated). However, 
the spatial momenta flowing along the propagators that are shown explicitly, and which lie within the dashed rectangles, are 
summed rather than integrated. 


up with an alternative to the very complicated and implicit definition given in Ref. Thus at this stage we only 
present the definition of T3, which is 


iF^ = 


iF 


1 ^ 1 
3 1 - iM 2 LiG 


iM2,LiF 


(3) 


Like F3, the quantities 1 /{ 2 ujL^), F, G, M.2,l are all matrices with indices k'fm';k£m. The hrst three of these 
quantities are known kinematic functions that depend on L (as well as i?, P, and m). They are defined in the 
appendix. The final quantity is the hnite-volume, two-to-two scattering amplitude, A 42 ,l- This in turn can be 
expressed (up to exponentially suppressed corrections) as 


iM2,L = 


1 

1 - I/CafP 


iK .2 , 


(4) 


where IC2 is a modified two-particle infinite-volume K-matrix introduced in Ref. [T]. The relation 0 is needed to 
show the equivalence of the result ([^ to the expressions for P3 given in Ref. [T] . 

Since M.2,l was barely used in Ref. [T], but plays an important role here, we discuss its properties in some detail. 
A useful, although imprecise, way of thinking about Ai2,L is as the two-particle scattering amplitude in finite volume. 
Recall that A 42 ,l das indices k't'm'] k£m, with k and k' short for k and k', respectively. For a two-particle quantity 
such as AI2.L, the momenta k and k' are those of the third, spectator particle, which is unscattered. Thus k = k', 
or, in matrix notation, cx 5 k'k- The same holds for /C2 and F. Thus the spectator momentum serves 

only to determine the energy-momentum flowing through the scattered pair, which is 


P2 = {E — uik, P — k), with ujk = y k'^ + rn? . 
The energy of the two scattered particles in their CM frame is therefore 

Elk = [P2] = [{E - ookf - (P - kf ], 

while the magnitude of their momenta in this frame is 

9^= [^^2*y4-m2]'/^ 


(5) 

( 6 ) 


(7) 
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The only unconstrained degrees of freedom are the incoming and outgoing directions of one of the particles in their 
CM frame, which we denote a* and d'*, respectively. 

We can now make precise the requirement mentioned above, namely that the theory should not have any two- 
particle resonances. The condition is that IC2 should not have poles in the range of allowed two particle energies. 
This requirement is most stringent for fc = 0 , since this maximizes the energy in the two particle subsystem. Thus IC2 
should have no poles in the range of two particle CM energies 0 < E2 q < E* — m. Such poles would lead to power 
law finite-volume effects when converting sums to integrals in diagrams like that shown in Fig. j^b), and such effects 
are not accounted for in the derivation of Ref. [T] . 

Returning to A^2,L) it is useful to have a diagrammatic definition in addition to the expression in terms of IC2 and 
E, Eq. Q. This is given in Fig.j^a), and defines fc’ Note that the spectator particle is not shown—it 

contributes an additional factor of Skk'- What appears in Fig.|^a) is the same sum of two-to-two connected Feynman 
diagrams as occurs in the infinite-volume two-particle scattering amplitude (with external propagators amputated 
and on shell) except that spatial components of loop momenta are summed rather than integrated. These sums can 
be replaced by integrals within the Bethe-Salpeter kernels (making only exponentially suppressed corrections), but 
not within the “boxed” loops. 

We now encounter for the first time an issue that will recur repeatedly in the following. This is that the external 
momenta in quantities such as A^2 ,l must be allowed to differ from finite-volume momenta (which we refer to as 
“not being in the finite-volume set”). By contrast, internal momenta are in the finite-volume set. This distinction 
is already present to some extent in Ref. [T], since when the non-spectator pair are put on shell their directions a* 
and d'* are maintained, but their magnitudes are changed to q^., so that they are no longer finite-volume momenta. 
Here we extend this by allowing d* and d'* to be arbitrary. This allows us to decompose the angular dependence in 
spherical harmonics and define the matrix form of A^2 ,l used in Eq. ([^: 

dTrYl,^,{d'*)M 2 ,L-k'e'm'-kimYe^{d*) = 5 k'kM 2 AElk,a'*,d*) ■ ( 8 ) 

Here and in the following there is an implicit sum over repeated indices. Note that M.2,l is not diagonal in its angular- 
momentum indices, unlike its infinite-volume counterpart, because rotation symmetry is broken by the restriction of 
internal momenta to the finite-volume set. 

In fact, we will need to consider more general external momenta when we take the infinite-volume limit in Sec. [V| 
This is because we take this limit holding the total momentum P, as well as the external momenta of individual 
particles such as fc, fixed. Thus, as L varies, P and k are not, in general, finite-volume momenta. This in turn 
implies that the total momentum flowing through the boxed loops in Fig. kla), which is P — k, is not, in general, 
in the finite-volume set. We define the extended M2,l by keeping the summed momentum in the set, i.e. of the 
form 27171/L, while the other loop momentum is not in the set. Since the Bethe-Salpeter kernels are symmetric under 
particle interchange, the choice of which momentum is summed is irrelevant. 

As noted above, the name “finite-volume scattering amplitude” is an imperfect moniker for A 42 ,l- For one thing, 
there are no in- or out-states in finite volume, and for another, Ai2,L is purely real (with no unitary cut). Nevertheless, 
we persist with the name because, as described in Sec. |Vj if we take the infinite-volume limit in an appropriate way, 
M2,l does go over to the infinite-volume scattering amplitude M2- 

In the following section we derive a result closely related to Eq. ([^ in which the left-hand side is replaced with a 
quantity that we call the finite-volume three-to-three scattering amplitude, and denote Ms^l- In direct analogy to 
M2 L, Ms L is defined as the sum of all on-shell, amputated six-point diagrams evaluated in finite-volume, as shown 
in Fig.l^b). As for M2,L^ the external momenta are not, in general, constrained by the finite volume, while some of 
the internal momenta are. The details are described in the following section. The quantities on which Ms^l depends 
are simply the degrees of freedom for two copies of three-particle phase space at fixed total energy and momentum, 
i.e. the dependence is Ms,L{k', d'*,k, a*). This is closely related to the dependence of M2,l{E2 a,, d'*,d*) except that 
now all three particles are involved in scattering, so the partitioning into “spectators” (fc and k') and “scattering 
pairs” (a'* and d*) is purely a choice of kinematic variables. In particular, it is no longer true that k' = fc, so that 
the CM energies and momenta of the scattering pair are different in the initial and final states. Nevertheless, one can 
separately boost to the scattering-pair CM frames, and define the directions of one of the pair as d* and a'* for initial 
and final states, respectively. This completely fixes the kinematic configuration. 

A key observation is that the diagrams contributing to Ms^l, shown in Fig. I^b), are closely related to those 
building up Cl, shown in Fig. Indeed, if we first discard the first diagram on the first line and the entire second 
line of diagrams in the latter figure, and then, for each remaining diagram, remove the tr and blobs, amputate and 
put on shell the six external propagators, we obtain exactly the diagrams for Ms^l- This observation allows us to 
make substantial use of the work of Ref. [I] . 
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{a) iAi2,L 
(b) iMs^L 
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+ 
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FIG. 4 . (a) The finite-volume two-to-two scattering amplitude, A4.2 ,l- We restrict attention to two-particle CM energies below 
four-particle threshold, < 4 m. As a result, only two-particle states can go on-shell and thus only loops shown explicitly 
in the figure have power-law finite-volume effects. The insertions between loops are two-to-two Bethe-Salpeter kernels, B2, 
evaluated in infinite-volume. As explained in Ref. pp, summing these diagrams leads to the result given in Eq. (b) The 
finite-volume three-to-three amplitude, Ada.n, defined as for M2,l but with six rather than four external legs. This contains 
both B2 and the three-to-three kernels, B^, which can be evaluated in infinite volume. S indicates symmetrization. In both 
(a) and (b), dashed lines attached to the kernels indicate amputated, on-shell propagators. Remaining notation as in Fig. 


As with A 42 ,l, calling a scattering amplitude for finite L is somewhat misleading, as it has no direct relation 

to scattering. However, as we discuss in Sec. does go over to the infinite-volume scattering amplitude Ms 

when L —>■ 00, as long as one takes the limit in the appropriate way. Thus, if we can derive a relation between 
and /Cdf,3; then taking the L —)■ 00 limit gives the desired relation between Ada and /Cdf,3- 


III. M3,l with B 2 KERNELS ONLY 


In this section and the next we determine the relation between and /Cdf,3. This section includes only the B2 

kernels in the skeleton expansions, while the next completes the job by adding kernels. This splitting mirrors the 
breakup of the derivation given in Ref. [T]. 

The contributions containing only B2 kernels are those of all but the first and last lines of Fig. Bb)- We find it 
convenient to work first with an unsymmetrized sum of these diagrams, which we labej^Adg'^j^’^^^^^^ (p, a'*\ k, a*). This 

quantity is defined by requiring that p and k are the momenta that are not scattered by the outermost B2 insertions. 
The superscript “m” indicates “unscattered”, and “(it, it)” means that this choice of coordinates is made in both 
initial and final states. This is the same notation as used in Ref. [ 1 ]. We decompose this function in harmonics in our 
standard way 






( 9 ) 


Note that is defined for any p and k, but in this equation we are restricting these to be finite-volume 

momenta. 

In Ref. [T], three-particle quantities similar to are needed only when P, p and k are finite-volume 

momenta. This implies that all three momenta in internal three-particle cuts in Fig. [^b) can simultaneously be in 
the finite-volume set. However, as noted in the previous section, when we take the infinite-volume limit we must 
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For brevity we have chosen the outgoing momenta to be p rather than k'. The definitions of a'* 
section. 


and a* 


are as for Ais^L iri the previous 
































FIG. 5 . Diagrammatic definition of Filled circles indicate K-matrices IC2, of which there are n. Solid lines are 

dressed propagators, the spatial momenta of which are summed (while time components are integrated). Dashed lines indicate 
amputated, on-shell propagators. Kinematic variables are defined as for 


consider general values of the external momenta. It is then not possible that all internal momenta are in the finite- 
volume set, and we must generalize the definition of This is straightforward for the diagrams in the top line 

of Fig.| 4 ];b): we simply take the two momenta that are being summed to be in the finite-volume set, implying that the 
third is not. Since B3 is symmetric, it does not matter which two momenta are chosen. In diagrams involving a boxed 
two-particle loop, we proceed as for Ai2,L, and choose (arbitrarily) one of the momenta to be in the finite-volume set. 
For the remaining figures, such as those in the third line of Fig.|^b), one has a non-trivial choice as to which of the 
loop momenta to sum. However, in all cases, the derivation of Ref. [T] provides a natural choice. The convention we 
use to define the extended is that the summed momenta always lie in the finite-volume set. We emphasize 

that the details of this choice are not relevant for our final result, since the distinction completely vanishes in the 
infinite-volume limit. 


To identify the desired egression for Ai-i,L we begin by introducing the relevant building blocks as presented in 
Ref. [T]. The first of these it[j 


iF 1 _ 1 iF 

2wL^ 1 - ilC2iF ~ 1 - iFi]C2 2ujL^ ' 


( 10 ) 


This quantity is related to M.2,l, given in Eq. (Ef,by 

'iA4.2,L = *^2 + {2u}L^)ilC2-A.iK,2 ■ 


( 11 ) 


Thus A contains the volume dependence of Ai2,L, entering through the factors of F. 

The second building block from Ref. [T] is a quantity defined in Fig. 5 This object is built from n > 2 

factors of JC2, with adjacent insertions scattering different pairs. It arises from secnons of diagrams where B2 kernels 
change at least once to a different pair of propagators. It contains power-law finite-volume effects, due to the presence 
of sums over the spatial momenta in loops. Summing over all possible numbers of IC2 factors leads to the quantity 


00 

y^{u,u) _ X '' y^{n,U,u) 

^3,L = / . ^3.L 

n—2 


( 12 ) 


Both K, 


and can be expressed as matrices in the same fashion as Ai 


3,L 


We can now state our result for the partial decomposition of A 4 


{u,u)[B2\ . 
3.L 


(u,u)IB2] 

3,L 


iAi 


{u,u)[B2\ 

3,L 


1 


1 - iK2iF 


00 


I 


n—0 


— iFilC2 


( 13 ) 


This decomposition makes explicit all finite-volume dependence arising from sum-integral differences (i.e. factors of F) 
on chains in which scattering occurs on the same pair of pr opa gators. It is a partial decomposition, however, because 
still depends on L. One can understand the result ( 13 ) qualitatively as follows: In order to have a connected 


scattering diagram, there must be at least one present. This can be “decorated” with chains of two-to-two 

scattering K-matrices and F factors on either the outside or between factors of . Summing the geometric series 

of two-to-two chains and then summing over any number of insertions leads to the result shown. 
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In Ref. ^ we denoted this quantity [,A], whereas here we drop the brackets. Note that l/{2uiL^) commutes with }C 2 and F, while JC 2 
and F do not commute. 
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If the external momenta are in the finite-volume set, then a derivation of Eq. (131 can be given by a straightforward 
application of the methods of Ref. [1]. We provide this derivation first, and then comment (toward the end of this 
section) on the generalization to arbitrary external momenta. 


The most direct path to a derivation of Eq. (13) is to start from Eq. (193) of Ref. [T], which we reproduce here: 




oo 


2 iF 

'3 2a;L3 ^ 

n=0 






(14) 


This is an expression for the finite-volume correlator including only B 2 kernels. Now, the main difference between 
Cl and is the presence of the “endcaps” a and cr^ in the former. When dressed by B 2 kernels these give the 

factors of and in Eq. (14), which are infinite-volume quantities. To obtain M.z,l we can discard these 

endcap factors by noting that the kinematic factor E, which contains the finite-volume dependence arising from a 
three-particle cut, amputates and puts on shell the adjacent matrices. Thus the factor between the outermost F 
cuts in Cl is precisely the amputated on-shell Ms^l that we are after, as long as the contribution contains two such 
cuts. There are no such cuts in the first two terms on the right-hand side of Eq. (14), since these are infinite-volume 


quantities, and so they do not contribute to M^^l- The third term also does not contribute since it contains only a 
single F cut. We conclude that we need only keep the final term, i.e. 


^/[B2 


A 


n—0 


(ic'xU 


^[52 


(15) 


If we drop the external factors of and and then further multiply by an inverse of iF/{2u;L^) on each side, 

we are left with all amputated six-point diagrams. This includes, however, disconnected diagrams (arising from the 
diagrams on the second line of Fig. |^. Since these are not contained in the definition of Ms^l they must be discarded. 
This is effected by summing from n = 1 rather than n = 0. In this way, and using the results in Eq. (10) for A, we 
obtain Eq. (13). 


In the remainder of this section we follow the method of analysis used in Sec. IVD of Ref. [T], where the result 


r D 1 

(14) for C\^ ^ is completely decomposed into finite- and infinite-volume quantities. Here we apply the method to the 


result (13) for In this case it turns out that we cannot simply borrow results from Ref. [T], but need to 

introduce two new quantities. Thus we provide a more detailed derivation. 

The starting point is to use the result for presented in Eq. (215) of Ref. [T], namely 


= iTiG[2ujL^]i]C2 + ( 1 iTiF ) |^i/Cdf.3 j ^ 



( 1 iTiF ) */Cdf,3 



^^T2u:L^ 


(16) 


where [Eqs. (212) and (213) of Ref. [T]] 


iT = 


1 


ziIC2. 


1 — iK,2iC 


‘'^df.3 
^^df.3 + *^df,3 


''^df,3 

' *^df,3 + *^df,3 


(17) 

(18) 


The result (16) separates finite-volume dependence—contained in F and C —from the infinite-volume quantities IC 2 


and /Cdf, 3 - All these quantities are matrices with our usual kim indices. In Eq. (18), /Cdf ,3 appears in various forms, 
with different superscripts, and it is convenient to collect these into a second level of matrices as shown. The differences 
between the forms of /Cdf ,3 are explained thoroughly in Ref. [T], and here we give only an overview. 

The simplest quantity to understand is /Cdf ’3 ^ • This is given by except for two changes: sums are replaced 

by integrals, and singular terms are subtracted (along the lines discussed in the introduction) so that it is divergence- 
free. The integrals use our non-standard pole prescription, and, in the case of multiple loops, must be ordered in the 
manner described in Ref. [1]. Since loop momenta are integrated, there are no subtleties defining for arbitrary 

external momenta. We stress that is not a physical quantity, because the symmetry between the three external 

particles is violated by our choice of coordinates. 
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To define the other quantities appearing in Eq. (181 it is necessary to change from the matrix form of 
which decomposes angular dependence into spherical harmonics, into a form with explicit dependence on the external 
momenta. This is achieved by 


IC^^ff{p,a'*-,k,an = ^7rY,l 


,[a J' 


(19) 


The other versions of /Cdf ,3 are then obtained by applying coordinate changes to this new form of We describe 

this explicitly for the example of R-ficall that, in the initial state, k is the momentum of the particle chosen as 

spectator (so that the first two-particle kernel B 2 scatters the other two particles) while a* is the direction of one of 
the other two particles in their CM frame. In the original frame the momentum of this second particle is called a —a 
momentum that is fully determined given k and a* (along with E and P). The third particle thus has momentum 
P — k — a in the original frame. A similar notation with k ^ p and o —> a' is used for the final state. As noted above, 
the first superscript of Adf ,3 is related to the choice of momentum variables for the final state, and the second to the 
choice for the initial state. Focusing on the initial state, quantities with superscripts 
choosing the spectator momentum to be a and P — k — a, respectively. For example. 


“s” and “s” are obtained by 


^df.3 


■,k,d*) = K!-^i^\p,a" 


,r) 


( 20 ) 


where k* is the direction of the particle with momentum k after it is boosted to the CM frame containing it and the 
particle with momentum P — k — a. Note that we could equally well have used the direction of the pa rticl e with 
momentum P — k — a after boosting in place of k* as the final argument on the right-hand side of Eq. (20), since 


^df ’3 ^ is symmetric under this interchange. Independently making analogous changes of variables in the final state we 


obtain all the versions of JCdi, 3 - We stress that this argument relies on being well defined for external momenta 

which are not in the finite-volume set, since, even if k and P are in this set, a is not. 

The final step to obtain the matrix forms of Adf ,3 appearing in Eq. (18) is to decompose the angular dependence 


in spherical harmonics. For example, the matrix form of 3 ^ is obtained from the left-hand side of (20) in exactly 
the same way as for [see Eq. (19)]. 


As is the case for all the other forms of /Cdf ,3 are separately unphysical. To obtain a physical quantity we 

must symmetrize over the initial and final assignments of momenta. In our matrix notation, this is achieved by 


( I 1 ) j l/Cdf,3 



( 21 ) 


Note that Eq. (16) does not contain the symmetrized version of Adf, 3 - A major part of the analysis in Ref. [T] involves 


r D 

rewriting in terms of the symmetric quantity, and the same is true in the following. 


To proceed, we substitute Eq. (16) into Eq. ( [T^ and organize the result in powers of /Cdf, 3 . For example, the 
contribution which is independent oDCdf.s is givenby substituting 


= iTiG[2ujL^]iK:2 + 0(/Cdf,3): 


( 22 ) 


[i.e. the first term in Eq. (16)] to find 


where 


iM 


{u,u)[B2] 
3,L 


= *I?i"’“'+0(/Cdf,3), 


(23) 


I LXJ 

{iTiG[2ujL^]iJC2A)''iTiG[2ujL^]ilC2 


1 - ilC2iF 


n—0 


1 

1 - iFilC2 


After some algebra using the definitions of A and T this simplifies to 


iV 


{u,u) _ 

L = 


1 


1 — iM.2,LiG 


T^*Al2,LiG*Al2,L[2wT^]. 


(24) 


(25) 
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The superscript on is a reminder that the external momentum indices are those of the particles not scattered 

by the outermost factors of A^ 2 ,l- 

Next we consider the term linear in /Cdf, 3 - To identify this it is necessary to substitute the Afdf, 3 -independent and 
(! 1 (/Cdf, 3 ) terms from Summing all contributions to the left and right of the single /Cdf ,3 insertion then gives 

[see also Eq. (224) of Ref. [T]] 




(«)’ 


(26) 


where 


( )^(i o) + - 


1 




- iM2,LiG 


iF 


— iM2,LlF{l 1) , 


0 ) “^ ( 1 y 2 wL3 1 - iM2,LiG 


-~iM2,L{‘^ojL^). 


(27) 

(28) 


These are the two new quantities that do not appear in Ref. [I]. The superscripts again indicate a particular choice 
of spectator momentum. 

To see the all orders pattern, it is sufficient to work out the C>(/Cdf 3 ) term. We find this to be 


iM 


{u,u)[B'2\ 


d( c 




‘3.L ^ \ '-L j 

where we have introduced two new “second-level” matrices 


iF-x 


2uj 


ilC. 


df,3 


7^ 


(tt) 


(^ ^ ) *^3 ( 11 ), 


IR. 

2u! 


^ \ iP.( 2 _1 

1 / 2w ^ 3 3 


(29) 

(30) 

(31) 


The factor between the /Cdf, 3 S in Eq. (29| is exactly the same as that in the corresponding analysis in Ref. [T] [see 


Eq. (230) of that work]. As in Ref. 


[[^ the term containing p is integrated, rather than summed, over spectator 


momentum (with sums over angular momentum indices remaining). More precisely, our notation means 




iAdf,3 = 




ilCdi,z{s,k) , 


(32) 


where p{k) is the phase-space factor defined in Eq. (A 6 ), and p and k on the right-hand side are finite-volume momenta 
which match the “outside” indices on the left hand side. The second level matrices involving /Cdf ,3 on the right hand 
side are the same as those defined in Eq. (18), except here we replace the spectator momentum indices with continuous 


momenta. We comment that the asymmetric form in which p{k) appears here arises from our convention for ordering 
the PV integrals, as is described in detail in Ref. [T]. 

Extending to all orders in /Cdf, 3 , and rearranging the sum, we deduce 




4“^ ) E 


n—0 


iFs 




[B2,p] 

df,3 


7^ 


(u)] 


(33) 


which is the analog of Eq. (234) of Ref. [T]. Here the p-dependent terms are summed into the new infinite-volume 
quantity 


^E(*^df,3 

n=0 \ / 




(34) 


We refer to these p-dependent terms as “decorations”. 
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The result ( |33[ ) has been derived assuming that the external spectator momenta and the total momentum take 
finite-volume values]^ We now argue that it holds for arbitrary choices of these momenta, as long as the finite-volume 
quantities F and G contained in and Af 2 ,l ars extended to arbitrary external momenta in the very straightforward 
manner described in the appendix. The derivation of Eq. (33) relies on repeated replacements of momentum sums over 
three particle cuts with momentum integrals plus the difference, One must ensure that replacements are done in 
an appropriate order (and that integrals are defined with the PV prescription), as discussed in detail in Ref. [1]. But 
what matters here are the following observations. First, the determination of which summands lead to a non-trivial 
sum-integral difference (i.e. the three particle cuts) is independent of the choice of external momenta. Second, the 
difference function F maintains its form (given in the appendix) for arbitrary external momenta. Third, the sums 
over spectator momentum indices, implicitly contained in the multiplication of matrices in Eq. (33) (and which are. 


by definition, over finite-volume momenta) correspond precisely to those loop momenta which should be summed over 
finite-volume momenta, according to the diagrammatic definition of given at the beginning of this section. 

In other words, the correct choice of which internal momenta to keep in the hnite-volume set is being made. And, 
finally, the kinematic factor G (which contains no momentum sums) extends to non-finite-volume momenta in a very 
straightforward way, as described in the appendix. Combining these observations, it is straightforward to see that the 
expression for given above remains valid for arbitrary external momenta, with one proviso. The proviso is 

that the external momenta, which propagate into the expression until there is a factor of G (either in or 

F^^) or a factor of are not in the finite-volume set. Thus the hrst matrix that appears on either end of the 

expression has different indices from matrices in the middle of the expression. In particular, the matrices on the end 
have one index which is an external momentum while the other is an (internal) hnite-volume momentum. Fortunately, 
these subtleties become irrelevant when one takes the inhnite-volume limit and sums become integrals. 


The hnal step in the construction of the contribution of from B 2 kernels is to symmetrize over initial and 


hnal momentum assignments. However, a complication arises because and F 

changes the assignment of momentum labels, the angular dependence now sweeps over the poles. Since this implies 
that the functions are not square integrable, one cannot decompose into spherical harmonics. To avoid this technical 
problem one must hrst change from spherical harmonic indices to an explicit dependence on angular variables, and 
only then symmetrize. For almost all choices of external momenta, one then avoids the singularity. 

To make the symmetrization operator explicit, we hrst convert to angular variables for in the standard 

way: 




(«) 


^ miLi ;v^ contain poles. When one 


A^K«)[B2] ( - ^ 


,(a'*)Al 


{u,u)[B2] 

3,L-,p£'m'-,kem^^rn{U ) . 


(35) 


Note that, as discussed above at length, we can treat all the external momenta as continuous variables, not constrained 
to lie in the hnite-volume set. We now dehne the quantities with u replaced by s and s in exactly the same way as 


discussed for /Cdf ,3 above—see Eq. 20 and surrounding text. We stress that this is simply implementing the interchange 
of the choice of spectator particle. Then we have 


x—u,s,s y—u,s,s 


(36) 

(37) 


The second line dehnes the action of the symmetrization operator S, and applies to any “(u, rt)” quantity that can be 
dehned for general external momenta. 

Applying S to both sides of Eq. (331 we hnd, after some algebra, the hnal result of this section. 


iA4 (p, a'* \k,a*) = iDl (p, a'* ',k,a*) 


S 


3 l-iM2,LiG 


iM2,LiF 


ilC 


[S 2 .P] 




-k 


iF 


3 2wL3 l-iM2,LiG 


iM2,l{‘^ujL^) 


■ (38) 


° The momenta of the scattered pair in the initial and final states are not in the finite-volume set, due to the on-shell condition, as 
discussed around Eq. The direction vectors for these momenta, a* and a'* ^ can be chosen arbitrarily. 

^ To see this most clearly, one should expand out the geometric series contained in ^ 3 , Ai 2 ,L in which case is a 

sum of terms each of which consists of products of either _F or G alternating with factors of /C 2 a-iid 
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Here is the symmetrized form of 




and 


^( 11 ) () (1 


(39) 


(40) 


equals the quantity of the same name appearing in Ref. [I] . 
Several comments are in order about the derivation of Eq. 


38[). When symmetrizing over the row vector 


(“) 


defined in Eq. (27l, the row vector (1 0) within it is converted into (1 1)13 Since the final form ( [3^ still contains 
an external symmetrization operator, we can replace (1 0) with (1 1) as long as we include a factor of 1/3. The row 
vector (1 1 ) can now be factored out and applied to the right. A similar discussion holds for with the symmetric 

column vector f ] factoring out at its left-hand end. Together with the fact that the matrix form of F^, Eq. (301, 


contains symmetric column and row vectors, this implies that all entries of the matrix in Eq. (33) are now 

sandwiched between symmetric vectors. This is why the symmetric quantity of Eq. (40) appears in our final 


result. One is then left with a geometric series of powers of which sums to the form shown in Eq. (js^. 

We stress that while symmetric under external particle exchange, violates this symmetry internally. This 

is because of the above-mentioned decorations involving the matrix p, which contains the asymmetric row vector 
(2/3,—1/3). Furthermore, our PV pole prescription breaks particle interchange symmetry. These facts make 
a difficult quantity to interpret. 


It is of course crucial that Eq. (38) contains the same g’^^ as appears in the analogous expression for given 

in Eq. (234) of Ref. [T]. We will see that this continues to hold when B 3 kernels are included. This means that the 
same quantity lCdi ,3 that appears in the quantization condition will appear in the relation to AI 3 . 


IV. M3 ,l with B2 AND Bs KERNELS 

In this section we complete the decomposition of M. 3 ^l by including all remaining diagrams, i.e. those with at least 
one B 3 kernel. This set is illustrated by the first and last lines of Fig. Bb)- Although the final result is straightforward, 
with the form of Eq. (38) being preserved aside from a new definition of /Cdf, 3 , the derivation of this result is quite 
involved. While we are able to partly reuse results from Ref. [T], significant additional work is required. As in the 
previous section, we derive the result first with finite-volume external spectator and total momenta, and then note 
that the result maintains its form if these momenta are allowed to be general. 

i ?3 kernels are symmetric functions of incoming and outgoing momenta that are smooth in the energy range of 
interest. They naturally divide diagrams contributing to into segments within which finite-volume dependence 

can arise, as illustrated in Fig.|^ These segments can be treated separately, and are of three types. The first connects 
the final state fields to the leftmost B^ kernel, and appears on the left-hand end of all diagrams considered in this 
section. The second type lies between two B^ kernels, and can appear any number of times (including zero) in the 
middle of the diagram. Finally, on the right-hand end there is always a segment connecting the rightmost B^ to the 
initial state. The second (“middle”) type of segment was analyzed in Ref. [T], while the other two are new to this 
work. 

We first recall the expression for a middle segment. To understand this we need the final result for finite-volume 
correlator including only B 2 kernels, . This is given by Eq. (239) of [T] : 


(^ 1 ® 2 ] _ pj[B2,p] _|_ ^ 


with 


Z = 


1 - 


iFs . 


(41) 


(42) 


Recall that our second-level matrix notation combines superscripts s and s into the lower entry; see, e.g., Eq. Jl8b. 
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(a) (5) (c) 


FIG. 6. Example of the decomposition of a contribution to M 3 ,l in terms of segments separated by insertions of the three-to- 
three kernel B 3 . As discussed in the text, the segments are of three types: (a) those connecting the final state to the leftmost 
kernel; (b) those which lie between two B 3 kernels; and (c) those connecting the rightmost kernel to the initial state. Segments 
of type (b) were analyzed in Ref. [I] , whereas those of type (a) and (c) are new to this work. 


This is obtained by performing manipulations analogous to those of the previous section on the starting form, Eq. (141. 

is the infinite-volume version of Cl with only B 2 kernels included, using the PV prescription, and including 
decorations due to p. It is defined in Sec. IVD of Ref. [1]. The second term in (411 contains the finite-volume 
dependence, which is collected into Z. All infinite-volume contributions between the cr at the left and the first finite- 
volume cut F are contained in while contains such contributions between the rightmost F and the 

ub Following Ref. [1], we define “decoration operators” Da' and Da such that [Eq. (241) of [1]] 



(43) 

A'[®2.p] = ^ ^ 

(44) 

A^B2,p] = d^B2,p]^^ 

(45) 


These are infinite-volume integral operators. The point of these definitions is that the decoration operators do not 
depend on the choice of functions at the ends. In these are a and but if both are replaced by B 3 then one 
obtains a middle segment of Thus the latter can be written 

... iBsM'^^^’PkBs , where , (46) 


which is Eq. (243) of [T]. 

We now turn to the end segments. We call the integral operator appearing in the left-hand segment L^'C[B 2 ,p]^ 
the corresponding operator for the right-hand segment r(“)[® 2 ’P 1 . Using this notation, the contribution to 
containing one and two RaS is 


No (u,u) superscripts are necessary for the middle segment due to the symmetry of B 3 . 

The result for r(“)[^ 2 >p 1 can be obtained by reflection from that for go we focus on the latter. Our starting 

point is a result that follows from Eq. (174) of [T] by the same argument given above for obtaining Eq. (13) from 
Eq. (pl). We find 


7^(ti)[B2,p]^t = 


1 


1 - i]C2iF 




n—0 


m—0 


(m,u) 

L 


3^ 


(48) 


Note that we are acting l(“)['® 2 ,p] qu cr^ rather than B^. This is just a convenience, as it allows us to reuse notation 
from Ref. [1]. The operator L^“) 1 ^ 2 ,p] jg the same for either choice. The new quantity is defined in Ref. [T], 

and represents diagrams in which is dressed with m factors of /C 2 j with the scattering pair switching between each 
factor. All loops except that closest to the cr^ involve sums over spatial momenta, and thus contain L dependence. 
This is illustrated in Fig. 17(a) of [1]. The lea ding term is = tjl (which will ultimately become R 3 in our case). 

The last term on the right-hand side in Eq. (48) is needed since the first term equals at leading order, but we 


require a"' /i in the unsymmetrized endcap. In the construction of this will correspond to factors of B^ being 

multiplied by 1/3 whenever they are the outermost insertions. This is the correct definition to recover a factor of B 3 
without the 1/3 from such diagrams within the symmetrized amplitude. Since B^ is itself symmetric, the effect of 
symmetrization on such diagrams is multiplication by 3, which cancels the 1/3 included here. 
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The form (48) is next rewritten using the results in Eqs. (188), (191) and (192) in [T], leading after some algebra to 


^ OO ^ 


(49) 


n=0 


Here and are infinite-volume quantities involving any number of 1 C 2 S connected to crb is the 

symmetrized version, which has already appeared in Eqs. (14) and ( [l^ . Both quantities are defined in Ref. HlH The 
result (49) is the analog of Eq. (13) above when there is a on the right-hand end rather than an initial state. 


Using the definitions of A and Eqs. ( |1()| and (13) respectively, we can rewrite our result for L^AIB 2 ,p] as 


3'^ 2a;L3' 


1 - iJC2iF 


aIB2 


(50) 


This allows us to reuse the work from the previous section, which led to the result (33) for Substituting 


the latter into Eq. (50), we encounter a new combination that can be rewritten as follows 


-T?U) 1 a\B2] _ 

2ccL3^ " 


iFs 


IR. 

2oj 


Al^^ 


where 


IR. 

2u! 


= ( 2/3 W . ^ ^ 
- ' -1/3 j 2a; ^ 


and 


^[52] = 


crt/3-f \ 

2crV3 + 24(")[^"l/ 

Here we have split up the symmetric quantity into a column matrix, such that 


^[B2l ^ ^ 1 1 W ^[Bs] I 


Introducing the definition [equivalent to Eq. (237) of [T]] 


^[B2,p] = -k 


I ' 


we find, after considerable algebra, that we can rewrite (50) as 




^[B2,p] 


n—0 


Finally, defining asymmetric decoration operators by 

/ £)(“)[B2.P] 

y p[B2,p\ _ ^(“)[B2,p] 




(51) 

(52) 


(53) 


(54) 


(55) 


(56) 


(57) 


(58) 


® In Ref. [T], denotes what we refer to here as Note that does not include an isolated factor of (j\ 

have at least one insertion of K. 2 - By contrast the symmetrized version, A^^^l does start with an isolated cr^ 


i.e. all terms in 

= at+C>(/C2)]. 
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we can write a stand-alone expression for the integral operator appearing in the left-hand segment 




iF:> 


n—0 


^(u)[B2,p] 


As for the middle segment, finite-volume dependence enters through the factors of F 3 . 

The corresponding expression for the right-hand segment is 


-R(“)[^2,p] = j' 25W[S2.P]^ jj[B2,p] _ J^{u)[B,,p] ) ^ J [ iF^ j [ 


n—0 


where we have defined the mirrored decoration operators 

fj ^ £){^)[B2,p\ £)[^2.p] _ jj{,u)[B2,p] ^ ^ ^'[B2.P] ^ ^ 

with 


7^ 




( 2l'[S2,p] ^ ^ ^ ^/[B2] ^ ^ ^^[B2] ^ 


2= 


where 


( A'l® 2 ] ) = ( cr/3 -h 2cr/3 + -h A'(") ) . 


(59) 


(60) 


(61) 


(62) 


(63) 


The quantities and are dehned in [T]. 

We now have all the pieces needed to construct the full unsymmetrized finite-volume scattering amplitude. The 
result is 


n—0 

Reorganizing the sums, and introducing the quantity 


(64) 


(65) 


where 



00 


E 

n=0 


^(n)[S2,p] 

£^[^2,p] _ ^^)[B2 ,p] 


iBsD 


[S2,P] 

c 


n 


iBs 


( 


r^[B2,p] 

B'a> ) 


jj[B2,p] _ ^(^)[^ 2 , P ] 


( 66 ) 


we find 


iM 


{u,u) 

3,L 


= iV 


'u,u) 


+ ( ^ 


» 



[B 2 &B 3 ] 

df,3 


00 


E 

n=0 


iFs 



[B 2 &B 3 ] 

df,3 



(67) 


As claimed above, this has exactly the same form as the result (331 from the previous section, but now the B 3 
contributions are contained within the modified /Cdf, 3 . As in the previous section, this result remains valid when the 
external momenta take arbitrary values instead of being restricted to the hnite-volume set. Indeed, the argument for 
this is unchanged, since the new features introduced by factors of B^ only impact the infinite-volume quantities in 
Eq. (67) for which the transition to arbitrary external momenta is not problematic. 

Symmetrizing proceeds as in the previous section, and leads to the final form for the finite-volume scattering 
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amplitudep^ 


iM3^L{p,a'*;k,a*) = WL(p,a'*;k,d*) 

1 


+ 5 

Here /Cdf, 3 , given by 


3 l-iM2LiG 


iM2,L'iF 


iFdt;. 


1 


l-iFsi/Cdf.a 


I ^ iF 1 

3 2wL3 l-iM2LiG 


iM2,l{‘^ujL^) 


*^df,3 ^ ( 1, 1 ) ( I) , 


• ( 68 ) 


(69) 


is equal to the quantity of the same name appearing in Eq. (251) of Ref. [I]. Thus we have achieved our objective— 
expressing l in terms of the same intermediate quantity /Cdf ,3 that appears in our quantization condition. 

The result (| 68 |) suggests an intuitive interpretation of /Cdf ,3 as an effective local three-particle interaction. In words, 
one can describe the second term on the right-hand side of Eq. ( 68 ) as having any number of two-particle interactions, 


jumping back and forth between scattering pairs, followed by a three-particle interaction (given by /Cdf.s), followed by 
more two-to-two scattering (T 3 ), possibly another three-particle interaction, etc., and ending with more two-to-two 
scattering. We note, however, that separating the local part of the three-particle scattering process involves some 
ambiguity, which here is manifested by the presence of the cut-off function H, whose detailed form is arbitrary. 


V. INFINITE-VOLUME LIMIT 


In this section we take the infinite-volume limit of Eq. ( 68 ) and thereby derive a relation between AI 3 , the standard 
three-to-three scattering amplitude, and /Cdf, 3 , the infinite-volume quantity appearing in the quantization condition. 

As a warm-up, we consider first the infinite-volume limit of the two-to-two finite-volume scattering amplitude, 
whose form we recall is 


■M2,L = 


1 

I+IC 2 F 


IC 2 


(70) 


with F defined in Eqs. (A3|-(A5|. M. 2 ,l, K ^2 and F are all matrices with both spectator momentum and angular 
momentum indices. The total energy and momentum is (E, P), while we call the incoming and outgoing spectator 


momenta k and p, respectively. As stressed repeatedly above, these momenta are not constrained to lie in the 
finite-volume set. We take L —>■ 00 holding these momenta fixed. 

The crucial question when taking the infinite-volume limit is how to treat sums over finite-volume momenta. We 
know that (suitably normalized) sums go over to integrals 



p 



(71) 


up to corrections which vanish as L —>■ 00 . If the summand is smooth the corrections vanish exponentially, while if 
there are singularities the corrections fall as powers of L. Keeping track of the latter was the main task in the analysis 
of Ref. [T]. Here, by contrast, we do not care how the corrections fall off, only that they vanish in the limit. However a 
subtlety arises for summands containing singularities, because the corresponding integrals are only well-defined with 
a pole prescription. For any finite L, the finite-volume momenta will avoid the singularities unless the energy E is 
equal to that of three non-interacting particles. But as L —>■ 00 , these singularities become arbitrarily dense and the 
limit is ill-defined. To take the limit we must first introduce a pole prescription for the sum. The natural choice is 
the ie prescription, since this is the choice used in defining scattering amplitudes in infinite volume. We also consider 
below what happens if we use the PV prescription. 


We note that there is a simple mnemonic that one can use to obtain this result from that for the full finite-volume correlator given 
in Eq. § This follows the procedure used to derive Eq. from the corresponding result for the finite-volume correlator, Eq. 

More precisely, identifying terms with at least two insertions of F, discarding the endcaps A' and A, multiplying each side by an 
inverse of 2 E/( 2 ajL^), discarding disconnected diagrams, and then symmetrizing, one reaches the final expression for the finite-volume 
three-to-three scattering amplitude given in Eq. ( |68[ >. This procedure cannot, however, be rigorously justified, for two reasons. The first 
issue is that various redefinitions of the infinite-volume quantities were required to reach Eq. § This obscures which diagrams are 
cont ribu ting to which quantities, and the only way we could resolve this uncertainty was by starting at an earlier stage of the derivation, 
Eq. | |14[ |, before such redefinitions have been made. The second issue is that the quantity reached by amputating A, JF as well as the 
outermost iFj has ill-defined exchange symmetry. It contains certain contributions which are symmetrized and others which are 
not. It turns out that symmetrizing the result in the manner described in the text removes this ambiguity, without overcounting or 
neglecting terms. 

















18 


This it prescription works for sums in a straightforward way. The procedure is as follows: (1) Replace all poles 
in sums by poles with nonzero ie; (2) Send L —> oo at fixed it —this limit is well defined with all sums becoming 
integrals; (3) Send e —>■ 0. In this limit the contribution of each diagram in the skeleton expansion for M 2 ,l [see 
Fig. I^a)] is converted into the corresponding contribution for M. 2 - Thus AI 2 .L —>■ M .2 when using this prescription, 
and similarly ^ Mz- 

Focusing first on A42,l, we note that it contains two types of finite-volume sums: those arising from matrix 
multiplications over the index k, and that contained in F’ru In this case (but not for M.z,l below) the former sums 
are trivial. This is because spectator momentum is conserved, so all matrices are proportional to Multiplying 

these out leada^to an overall 5kp- 

The non-trivim sum is that contained in F. Naively one might say that, since F contains a sum-integral difference 
[see Eq. (A5|], F ^ 0 when L ^ 00 . After all, F is the source of hnite-volume corrections. However, as already 
noted, this conclusion depends on the pole prescription chosen to define the sum before the limit is evaluated. If we 
choose the PV prescription for the sum, then, indeed, F vanishes in the limit 


lim 

L—¥00 


E = 0, 


(72) 


PV 


Here we denote the choice of pole prescription by the s ubsc ript. On the other hand, if we dehne the limit by regulating 
the sum with an it prescription, then it is of Eq. (A5| that vanishes, so that 


lim 

L—>-oo 


F — ^k'k — P 


(73) 


where is the phase space matrix given in Eq. (A6). We stress that, in the following, p without an argument 


is a matrix with full k£m indices while p{k) has only im indices. 

Using these results we can immediately determine the infinite-volume limit of M 2 ,l- If we use the PV prescription 
for the limit of sums, then 


lim 

L—¥CCi 


AI2.L — ^2 


(74) 


PV 


Above threshold, PV is equivalent to the principal value prescription. Thus Eq. (74) reproduces the well known result 


that evaluating scattering diagrams using the PV prescription in loops (i.e. dropping the imaginary part) leads to 
the K-matrix. This is not, however, the infinite-volume limit that we desire. If we instead use the it prescription we 
obtain 


lim 

L—^OO 


■M2,L = M.2 = 


1 


\FK2P 


-K 2 


(75) 


Above threshold, this is indeed the standard relation between M .2 and 1C2, with p adding in the cuts needed for 
unitarity. This provides a check that we are taking the infinite-volume limit correctly. Below threshold (i.e. with 
the spectator momentum chosen so that the scattering pair are below threshold), this result provides an alternative 
expression for /C 2 in terms of A42 analytically continued below threshold: 

^t=M-^^-p. (76) 

Note that, although /C 2 and p depend on the cut-off function H{k), this dependence cancels in Al 2 - 


We now apply the same limiting procedure to our result for M.z,l^ Eq. (68). To do so we need some new notation. 
Since p and k become continuous variables, while angular-momentum indices remain unchanged, we adopt a hybrid 
notation for /Cdf ,3 and e.g. 


/C 


3,df-,pi'm' ikim 


= /c 




(p, k). 


(77) 


The matrix multiplications over the angular momentum indices im go over unchanged to the infinite-volume limit. 

There is a minor inconsistency between our matrix notation—designed for summing over finite-volume momenta—and the fact that 
we have generalized to arbitrary spectator momenta. In particular, since the external momenta do not appear in the finite-volume 
set, what does the mean? Working back through the derivation, we find that it is simply a mnemonic for spectator-momentum 

conservation, with no sums over intermediate spectator momenta actually appearing. The same comment applies below for all instances 

of this Kronecker delta. _ _ 

This is because the p term in Eq. | |A4[ | is exactly what is required to turn the integral with the ie prescription in Eq. into an 

integral with the PV prescription. See Ref. ^ for further discussion. 
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We stress here that /Ca^df was already defined for arbitrary p and k above so that Eq. (77) is only a relabeling. 
Similarly, we make the continuous spectator momentum an argument of AI 2 , and drop the Kronecker delta from its 
definition: 


lim 

L—^oo 


^klm — rn'\irn{jP^ ■ 


(78) 


This is similar to the notation for p in Eq. (73). Finally, for AI 3 and its divergent part T), which, after symmetrization, 
are functions of angular variables rather than spherical harmonics [as discussed in Sec. Ill around Eq. (36)], we use 


lim 

L—>-oo 


M 3 ^L{p,a'*-,k,a*) = M^{p,a'*]k,a *), 


(79) 


and similarly for T). 

Our task is now to evaluate the limit in 


iAi 3 {p,a'*]k,a*)= lim 

L—>-oo 


1 


si iV 


(u^u) 


3 l-iM2LiG 


iM2,LiF 


iAfdf,: 


1 




I ^ iF 1 

3 2wL3 l-iM2LiG 


zM 2 ,L( 2 wi^) 


(80) 


We can interchange the symmetrization and the taking of L —>■ 00 because the changes of variables involve d in 
symmetrization are independent of L. We thus begin by determining the L —)■ 00 limit of given in Eq. (25). 

To do so we replace M. 2 ,l with AI 2 , and introduce an infinite-volume version of G 




^T:Ytm.'(k*)H{p)H{k)Yl^{p*) fp 


2u!kp{E - uJk - ujp - ujkp + ie-) \ql 


(81) 


This differs from the matrix form of G, given in Eq. (A2), simply by the removal of a factor of l/{2ujkL^) ■ This factor 
either cancels the explicit {2ujL^) in or converts sums over internal spectator momentum indices into integrals. 

For example the first two terms in the geometric series obtained by expanding have the following limits: 


lim 

L—>oo 


iM2iGiM2{2ujL^) = iM2{p)tG°°{p,k)iM2{k) 


lim 

L—>-oo 


iM2iGiM2iGiM2{2ujL^) = / —iM2WG°°{p,^iM2WG°°{s,k)iM2{k), 

./c 


(82) 

(83) 


where angular momentum indices are implicit. Here we are using f^ = f d^s/(27r)^. Repeating for higher-order terms 
we find that the infinite-volume form of 


lim 
L—>^00 


^L\pl'm' \klm 




{p, k) , 


(84) 


satisfies the integral equation 

iV^'^’^\p,k) = iM2{p)iG°^{p,k)iM2{k)+ [ ^iM2{p)iG'^{p,s)iV^'^’^\s,k). (85) 

Js 2 ws 


Symmetrizing following Eq. (37) leads to 


I?(p,a'*;fc,a*) =5{l?(“'“)(p,fc)} . ( 86 ) 

This quantity contains all the physical divergences in AI 3 , i.e. divergences that occur above threshold for physical 
external momenta. The difference 

Mdf, 3 iP,d'*;k,a*) = M 3 {p,d'*;k,d*) - 'D{p,d'*;k,d*) (87) 

is free of such divergences, and can, if desired, be expanded in spherical harmonics. 
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To determine the L ^ oo limit of the remainder of we need to evaluate 


lim 

L—>oo 


1 


1 


3 l-iM 2 LiG 


iM2,LiF 


iK-dt- 


1 


1 


iF 


1 


1 — iF^ilCdf^s 

where F^ is defined in Eq. ([^. We first take the infinite-volume limit of the central part of (881, i.e. 

1 


lim 

L—^oo 


iFd{,3 


1 - iFsilCdfR 


3 2a;L3 1 - iM 2 ,LiG 

t of the central pa 

= ^'T£'rn';£m{P: • 




Rewriting F 3 as 


iFd = 


iF 

2wL3 


- p£'m;k£m 




we find that T is given by the solution to the integral equation 


iT{p,k) = ilCdi,3{p,k) + I / i/Cdf,3(p, fc), 


where 


£(“.“)= ( - + iM2ip)ip{p)) {2TrfS^{p- k) + iVG’-^\p,k) 


.ipjk) 

‘2uJk 


( 88 ) 


(89) 


(90) 


(91) 


(92) 


and one must additionally enforce that T is symmetric under particle interchange symmetry on its left-hand argu¬ 
ment]^ The factors of p{k) arise from the infinite-volume limit of f 

What remains to be considered are the combinations at the left- and right-hand ends of Eq. (88). These become 
integral operators in the infinite-volume limit. Putting everything together we find 


iM 3 {p, a '*; fc, a*) = iV{p, o'*; fc, a*) + 5 I ^ ^ (p, s)ir{s, f )7^(“’“) (r, ) | , 

where is the reflection of 

7^(“’“^(p,fc) = Q F ip{p)iM2(:p) \ (27r)^a^(p- k) + (p, fc). 


(93) 


(94) 


Eq. (93) is our final result. It shows that M 3 can be obtained from /Cdf ,3 by solving two integral equations (first 
for which requires A42, and then for T), applying two integral operators to T, symmetrizing, and then adding 

in the divergent part T). We note that all the integrals appearing above have a finite range because all integrands 
include the cutoff function H. This may be helpful when making numerical approximations. 


VI. EXPRESSING Vdf 3 IN TERMS OF Ma 


In this brief section we show how the result (93) can be inverted so as to obtain Afdf.a given M 3 (and M 2 )- This 


will be useful if one has a model for M 3 and wants to determine the prediction for /Cdf, 3 i which can then be inserted 
into the quantization condition of Ref. [T] in order to predict the finite-volume spectrum. 

We begin by rewriting Eq. (93) in terms of symmetric quantities: 


tMdiAP,a'*-,k,a*) = j[i27T)^S^p- s) 47 rS^{a'* - b'*) + Ac{p,a'*-,s,b'*)} 

xiT{s,y*;f,b*) \^{2Tr)^S^{r- k)4TrS^{b* - a*) + An{r,b*;k,a*)'^ . (95) 


One way to do this would be to symmetrize the matrix appearing between /Cdf 3 and T. 

One would expect the final result to involve integrals over Lorentz-invariant three-body phase space, and indeed this is the case. The 
implicit sum over angular momentum indices can be rewritten as an angular integral in the scattering pair’s CM frame. This integral, 
combined with that over the spectator momentum, gives three-body phase space (extended to include subthreshold particles) as long 
as one includes the Jacobian factor p(fc)/2ajfc. 
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Here the angular averages and corresponding delta-functions are defined by 

Iri 


dflr, and 4TrS^{b — d)f{b) = f{d). 


(96) 


These angular averages arise from reexpressing the sums over repeated £m indices as angular integrals. We have 
also changed from spherical harmonic indices to explicit angular dependence for the quantity T, using our standard 
normalization [see, e.g., Eq. (35)]. Finally, the quantity At^ is defined by 


{27r)^S^{p-k)4Tr6^{a'*-d*) + Anip,a'*-,k,d*) = ls{TZ^^’^\p,k)} , 


(97) 


with Ac defined analogously. Note that since 77.^“’“) and £(“>“) are singular quantities (due to the delta functions 
and the poles in they must be symmetrized in the manner described by Eq. (37). This is why we must change 

to angular variables. 

To obtain Eq. (95) from Eq. (93) we have, apart from the change to angular variables, used the fact that T is 
symmetric under independent particle exchange on both its initial and final arguments. This allows us to symmetrize 
and on both the outer and inner arguments, which is what is done by the operation of 5. To avoid over 

counting we must include the factor of 1/3. We note that, although and At^ are complicated quantities, one can 
in principle construct them given A^ 2 - 

We now invert Eq. (95) in stages. The first step is to construct A^df .3 from AI 3 using Eq. (87). Next we invert the 
factors on the left and right sides of T in (95). This is done using kernels Ic and In that solve the integral equations 

Ic(p, d '*; s, b*)Ac (s, b*;k,d*), 


Icip,d'*]k,d*) = {2tt)^S{p— k)4'KS'^{d'* — d*) — 
lTi{p,d'*;k,d*) = {2tt)^S{p— fc) 47 r 5 ^(d^* — a*) — 
One then finds that 


s Jb* 


A-nip, a*; s,b*)In{r,b*; k, d*). 


s Jb’*’ 


{d'*)lT{p, k )pm-,lmytm{d*) = 


[ Ic{p, d'* ; s, 6 '*)*Mdf ,3 (s, b '*; f, b*)In (r, b*-k,d*), 


(98) 

(99) 

( 100 ) 


where we have used the fact that T is non-singular to keep it expressed with spherical harmonic indices. 
Finally, we reconstruct TCdf.s from T by rewriting Eq. (91) as 

iK^dfMP^k) = iT{p,k) - I / z/Cdf, 3 (p,s)*|^(s, f)jT(r, fc), 


( 101 ) 


and then solving for /Cdf, 3 - This must be done with the added condition that /Cdf ,3 is symmetric under particle 
interchange symmetry. 


VII. SIMPLIFYING CASES 


In this section we describe two approximations in which the relation between Ads and /Cdf ,3 simplifies substantially. 
Such simplifications are likely to be the first step towards using this formalism in practice. Both approximations were 
discussed in Ref. [T]. 

We begin by considering the result in the approximation that only the s-wave components of AI 2 and /Cdf ,3 are 
nonzero—which we call the s-wave approximation. In this case the quantities /Cdf, 3 , T, C, IZ and Ai 2 all are 

non-zero only for £ = to = 0, so the spherical harmonic indices can be dropped. This implies that Eqs. ( |85[ ), (91), 
(92) and (94) retain their forms, except that now all quantities are no longer matrices. In addition, becomes the 


function 


G“(p,fc) = 


H{p)H{k) 


‘2dJkp{Id OJk djp UJkp iY) 


The final equation (93) is replaced by 


iAis{p,d'*\k,d*) = S^^V^’'^\p,k) + J J £*^“’“^(p, s )zT(s, r (r, fc) 


( 102 ) 


(103) 
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where symmetrization is defined as before [Eq. (37)] except that now the quantity being symmetrized has no depen¬ 
dence on d'* and d*. Nevertheless, such dependence is introduced by the symmetrization, and thus Ais maintains its 
dependence on the full set of kinematic variables (albeit in a simplified form). 

Next we consider the more restrictive approximation—referred to in Ref. [T] as the isotropic approximation—in 


which the s-wave approximation is extended by assuming that lCdf, 3 iPy depends only on the total three-particle CM 
energy E*, 




(104) 


In this case, it is shown in Ref. [1] that can be obtained from the spectrum using a simple algebraic equation 

[Eq. (38) of P]. The integral equations which determine AI 3 are all at fixed E and P and thus at fixed E*. Thus 
simply a constant, which simplifies some of the results. 

We find that is the same as in the s-wave approximation, while C and TZ simplify to 


7^(fc) 


Q^+iM2ip)ip{p)^ -t 

+ ip{k )iM 2 {k )) + ^ (C 


(105) 

(106) 


Introdncing the new quantity 




■ipjs) 


iM2{s)ip{s) 



'ipjs) 

- 




f) 


ip(^ 

2uJr - ’ 


(107) 


the final result is 


iAid{p, d'*; fc, a*) = S 


iV 


{u,u) 

3 


{p, k) + 


1 



(108) 


Thus, in this approximation, the only integral equation that needs to be solved is that for Apart from that, 

the result for AI 3 is obtained simply by doing integrals, all of which are over a finite range due to the cutoff function 


H. 


VIII. CONCLUSIONS AND OUTLOOK 

In this work we have derived a relation between the non-standard three-particle K-matrix, /Cdf, 3 , and the standard 
three-to-three scattering amplitude, This completes the formalism relating the finite-volume three-particle 

spectrum to scattering observables and removes the central drawback of the formalism presented in Ref. [T]. 

Euture work will be dedicated to providing non-trivial checks of the formalism completed here, as well as assessing 
its utility. Concerning the former aim, we are currently completing a study of the three-particle quantization condition 
for P = 0 and for energy close to three-particle threshold, E « 3m. Here we can compare our result to those obtained 
using non-relativistic quantum mechanics mm. and using perturbation theory in a theory [7]. As we describe 
in our upcoming work, we find nontrivial agreement between our general result in this limit and the results of the 
alternative analyses. 

Concerning the utility of this framework, one important focus will be on solving the integral equations relating 
Adi ,3 and AI 3 . A possible direction that one might pursue follows from the observation that, with total energy and 
momentum fixed, the phase space for any number of particles is compact. This means that smooth functions which 
depend on this phase space, such as M 2 in the case of two-particles and Adi ,3 in the case of three, must be expressible 
as an infinite series of generalized harmonics. In the two-particle case this is the standard decomposition in spherical 
harmonics, leading to the parametrization of M 2 in terms of an energy-dependent phase shift for each partial wave. 
Identifying an analogous decomposition for Adi ,3 would be a first step in providing a well-motivated parametrization 
to extract this quantity from the spectrum and would likely also simplify the procedure for constructing Ms- 

Einally, we stress that the most important function of this work is as a stepping stone to the formalism required 
to describe more complicated and more physically interesting systems. We thus aim, in future work, to remove the 
restrictions listed in the introduction and provide a formal result for all possible three-hadron states. Eor example, to 
accommodate nucleon channels we must incorporate non-zero intrinsic spin into this analysis. This has been studied 
thoroughly in the two-particle sector in Refs. [IIHS]. We expect the methods developed there will facilitate the 
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generalization to three-particles. A further necessary extension is to accommodate two-to-three transitions, as well as 
transitions involving multiple two- and three-particle channels, possibly containing non-identical and non-degenerate 
particles. Although this will introduce many complicating details to the analysis, we expected that much of the 
technology developed here can be reapplied. Finally, we hope to extend this work to apply to systems with two-body 
bound states and resonances. This will require accommodating new finite-volume effects arising from bound-state 
and resonance poles. 
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Appendix A: Definitions 

We collect here definitions from Ref. [1] that are not given in the main text. The matrices appearing in F 3 , Eq. ([^ 


are 


[2ujL^ 


I k'£'m' :kim 




Gpt 


m':k£m — 




€ 


^k ^kp') 

^k'£'m'\klm — ^k'k^i'm'i 


1 


2ujkL^ ’ 


Flfm'-,im{k) = ^ 




d^p 

(27r)3 


inYi,rn'{p*)yi*rniP*)Hik)H{p)H{hp ) f p* \ 
2ojp2uJkp(^E ojf^ iOp cc/cp T 


£+i' 


<ll) 


(Al) 

(A2) 

(A3) 

(A4) 

(AS) 


where the sum over p in F runs over all finite-volume momenta. Here bkp = P — k — p is the momentum of the third 
particle, and ojkp = \Jb\p + 'm? is its on shell energy. The momentum p * is the result of boosting p to the frame in 
which the non-spectator pair [which has four-momentum P 2 , see Eq. ([^] is at rest. If all three particles are on shell, 
then jp *1 = p* = < 7 ^, where is given in Eq. (j^. The quantities k* and q* are obtained similarly by interchanging 
the roles of k and p. The matrix p is a phase-space factor for the non-spectator pair, defined by 


P£'m’\£m{k') — d g’tdm’mH {k^p{P2) 


m) = 


j—i^/P^J4YFrrP ( 2 m)^ < P|, 


167 rv/if \ I VF|/ 4 -™ 2 | 0 < P| < ( 2 m )2 . 

Here, and in the definitions of G and F, H provides a smooth ultraviolet cut-off. We require 


H{k) = 


0, P|<0; 

1 , ( 2 m) 2 <P| 


(A 6 ) 

(A7) 

(AS) 


where the first condition removes unphysical boosts and the second ensures that the cut-off does not change the 
contributions from on-shell intermediate states. In the intermediate region, 0 < P| < (2m)^, H{k) interpolates 
smoothly between 0 and 1. An example of a function which satisfies all requirements is 


H{k) = JiPi/[4m ^]), 


with 


' 0 , 

J{x)= { exp (-y exp 

. 1 , 


1 

l — x 


X <0; 

0 < a; < I; 
I < a;. 


(A9) 


(AIO) 
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The definitions in Eqs. (A1)-(A7) are used in Ref. [T] solely with the spectator momenta {k, k' and ^ as well as the 
total momentum (P) all being finite-volume momenta. As noted repeatedly in the main text, in this work we need 
to consider the extension to the case that these momenta are not in the finite-volume set. This presents no problems 
for any of the above definitions—all extend smoothly to general values of k, fc', p and P. 
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